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I ^' Abstract 

^ ■ 

Q^' The violation of J. Bell's inequality with two entangled and spatially separated quantum two- 

pH , level systems (TLS) is often considered as the most prominent demonstration that nature does 

c^ : 

^ I not obey "local realism". Under different but related assumptions of "macrorealism", plausible for 

macroscopic systems, Leggett and Garg derived a similar inequality for a single degree of freedom 

^ ■ undergoing coherent oscillations and being measured at successive times. Such a "Bell's inequality 

^n : 

ff^ , in time", which should be violated by a quantum TLS, is tested here. In this work, the TLS 

Cn I is a superconducting quantum circuit whose Rabi oscillations are continuously driven while it is 

wS ■ 

^^ ■ continuously and weakly measured. The time correlations present at the detector output agree with 



o 



quantum-mechanical predictions and violate the inequality by 5 standard deviations. 



^ ! PACS numbers: 74.50,03.65,82.25 



Introduction 

The violation of J. Bell's inequality l|, |2| is the most prominent example of a situation 
where the predictions of quantum mechanics are inconipatible with a large class of classical 
theories. In the early 1980s, Aspect and coworkers [3] brought an experimental proof of 
this violation using pairs of spatially-separated polarization-entangled photons. By demon- 
strating an excess of correlations between the polarizations measured on the two photons of 
a pair, they ruled out descriptions of nature satisfying the very general conditions known 
as local realism. This striking finding also contributed to transform the so-called quantum 
weirdness into a useful resource for information processing. Shortly after, quantum cryptog- 
raphy protocols and quantum algorithms exploiting entanglement were indeed proposed |4|. 
Following a reasoning similar to that of Bell, Leggett and Garg derived in 1985 an inequality 
that can be seen as a "Bell's inequality in time", which applies to any single macroscopic 
system measured at successive times |5| and fullfiling the assumptions of macrorealism: (Al) 
the system is always in one of its macroscopically distinguishable states, and (A2) this state 
can be measured in a non invasive way, i.e. without perturbing the subsequent dynamics of 
the system. Quantum mechanics however contradicts both assumptions, which can lead to 
an excess of correlations between subsequent measurements and to a violation of this inequal- 
ity. Ruskov and coworkers [6| then adapted the inequality to the situation where a two-level 
systems (TLS) is continuously and weakly monitored during its coherent oscillations. Us- 
ing such a weak monitoring, we report here an experimental test of a Bell's inequality in 
time (see also the recent works 7|, |8|), yielding results in excellent agreement with simple 
quantum-mechanical predictions and in contradiction with a large class of macrorealistic 
models. 

Bell's inequalities in space and in time 

We start by briefly recalling the experimental protocol of the usual CHSH test J2|] of 
Bell's inequalities (see Fig. [1^). It consists in identically preparing many times a pair of 
quantum TLS in a maximally entangled state such as l^J') = (|ti) ~ lit))/"\/2- Each 
member of the pair is then distributed to two observers A and B, who perform projective 
measurements of the TLS spin 0"^ ' = ±1 along one of two directions Oj (i = 1, 2) for A and 
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Figure 1: Comparison between two thought experiments which test the usual CHSH Beh's inequahty 
and the Beh's inequahty in time. (A) CHSH inequahty: two maximahy entangled spins aA and as 
are sent to two spatially separated observers A and B. Each of the observers measures with pick-up 
coils his spin along one of two possible directions (oi and 02 for A, and 61 and 62 for B ); the 
four directions make angles 6 as depicted. By repeating this experiment on a statistical ensemble, 
a linear combination S of the four possible correlators between measurements on a spin pair is 
computed. Local realism requires —2 < S < 2, while quantum mechanics predicts S = 2 -v/2 for 
9 = 45°. (B) Bell's inequality in time with weak measurement: a single spin a undergoing coherent 
oscillations at frequency w/j is continuously measured with a pick-up coil coupled to it so weakly 
that the time for a complete projective measurement would be much longer than the period of 
oscillations Tr = 2-11/1011. From the noisy time trace recorded in the steady state, one computes 
a linear combination fic of the three time- averaged- correlators between the readout outcomes at 
three times separated by r. Macrorealism requires /lg ^ 1 foi' any r, while quantum mechanics 
predicts f^c = 1.5 at r = Tr/6. 



bi for B, with these directions forming angles (ai,6i) = (61,02) = (02,^2) = ^ as shown in 
Fig. [1^. The two observers then combine all their measurements to compute the Bell sum 
E{e) = -Ku + K12 - K22 - K21 of the correlators Kij{e) = {afaf). The Bell's theorem, 
based on a simple statistical argument, states that according to all local realistic theories 

-2<S(^)<2. (1) 

However, standard quantum mechanics predicts that this inequality is violated, with a max- 
imum violation T,{9 = 7r/4) = 2v2. Many experimental tests, and in particular those 
performed by A. Aspect [3;] have verified this violation jlll. Il2l|. 

While quantum entanglement between two spatially separated TLS is at the heart of the 
previous violation, Leggett and Garg proposed a similar inequality |5i] holding for a single 
degree of freedom — 1 < z{t) < 1 fulfilling the assumptions of macrorealism (z(t) defined at 
any time, and measurable with no perturbation). Using a simple arithmetic argument a la 
Bell, they showed that 

z{to)z{h) + z{h)z{t2) - z{to)z{t2) < 1 (2) 

for all {tj}. Consequently, an observer measuring z on many identical systems, either at to 
and ti = to + T, or at to and t2 = to + 2r, or at ti and t2 should find ensemble-averaged 
correlators Kij(to,T) = {z{ti)z(tj)) (for i,j = 0,1,2, with i < j) satisfying the Leggett- 
Garg's inequality: 

hcito, r) = Koi + Ku - K02 < 1. (3) 

Quantum mechanics on the other hand predicts that, applied to the case of a quantum 
TLS undergoing coherent oscillations at frequency Uji , this inequality is violated for well- 
chosen values of r, with maximum violation /lg(^0;''" = '^/'^^r) = 1-5 independent of to. 
Here, the delay r between successive measurements plays the role of the angle 6 between the 
measurement directions in the Bell's inequality ([1]), justifying the nickname "Bell's inequality 
in time". The excess of correlations predicted by quantum mechanics, compared to the 
macrorealistic case, can be interpreted as resulting from the projection of the TLS state on 
a az eigenstate induced by the first measurement. 

As shown in [6], the very same conclusions also hold if the TLS undergoing coherent 
oscillations is continuously and weakly monitored along cr^ (see Fig. IB) instead of being 



projectively measured at well-defined times. The detector now delivers an output signal 
V{t) = {6V/2)z{t) +C,{t) proportional to z(t) with some additional noise ^(t). Macrorealism 
implies that the dynamics of the system at time t + r is fully uncorrelated with the detector 
noise at time t so that {C,{t)z{t + t))j = 0. The detector's output correlation function 
K{t) = {V{t)V{t + T))t/{SV/2Y is then simply equal to {z{t)z{t + r))^. By averaging 
inequality ([2]) over to in the steady-state, the Bell's inequality in time ([3]) becomes 

fLG{r) = 2K{t) - K{2r) < 1, (4) 

and should be violated by a quantum TLS in the very same way as discussed above. Here 
the violation is however not due to a strong projection of the TLS wavefunction induced 
by measurements at well-defined times of its evolution, but rather to the continuous partial 
projection caused by the measurement during the TLS coherent evolution, which reinforces 
correlations between the detector output at successive times. 

Experimental setup 

Our experimental setup (see Fig. [2JA. and supplementary informatior T^ for probing in- 
equality (jlj) closely implements the proposal discussed above, making use of the possibilities 
offered by the so-called circuit quantum electrodynamics (circuit-QED) architecture 13|, [lj| 



where a superconducting artificial TLS is coupled to a superconducting coplanar waveguide 
resonator. The TLS consists here of the two lowest energy states g and e of a modified 
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Cooper-pair box of the transmon type 17|, [iSJ. These two states can be regarded as "niacro- 
scopically distinguishable" because the dipole moment of the g — e transition is of the order 
of 10^ atomic units. On the other hand, the only degree of freedom of this system is the 
phase difference between the superconducting order parameters on both sides of the Joseph- 
son junction forming the Cooper-pair box, conjugate to the number of Cooper pairs passed 
through the junction; this phase is a collective variable whose degree of macroscopicity is 
still under debate (iSl. Il6|. 

The TLS transition frequency is Uge/^ix = 5.304 GHz, below the resonance frequency 
Uc/2'K = 5.796 GHz of the resonator to which it is capacitively coupled for its measurement. 
Two microwave sources Vd and Vm drive and measure the TLS at frequencies Uge and Uc, 
respectively. In order to continuously monitor the induced Rabi oscillations up to a few 
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Figure 2: (A) Experimental implementation of the thought experiment in Fig. [TJ3) with a quantum 
electrical circuit. The spin (or TLS) is a Cooper pair box of the transmon type (magenta) capaci- 
tively coupled to a microwave resonator (cavity sketched as a green coaxial cable) . Two microwave 
sources V^ and Vm are used to drive and measure the transmon at oOge and cjc, respectively. The 
reflected microwave at uj^ is routed through a circulator to a cryogenic amplifier followed by an 
I-Q demodulator. The time dependent phase ip{t) measured by the demodulator carries the in- 
formation about the TLS state. (B) Measurement-induced dephasing of the TLS as a function of 
the measurement strength, i.e. the amplitude of V^ or the mean photon number n in the cavity. 
Measured ensemble-averaged Rabi oscillations (c/9(At)) in presence of n = 0, 1, 2, 5, 10 and 20 
photons (from blue to red curve), showing the transition from weak to strong measurement. Each 
curve is fitted to a solution of Bloch equations (thin black lines), in quantitative agreement with 
expected measurement-induced dephasing rates (see online supplementary informatior fP]) . 

tens of MHz, we implement a resonator bandwidth of kJI-r = 30.3 ± 0.8 MHz (quality 
factor 191 ± 5) by designing the appropriate resonator input capacitance [lj|. With these 
parameters, the TLS is sufficiently detuned from the resonator for their interaction to be 
well described by the so-called dispersive Hamiltonian H = hxndz ll3|, with n the photon 
number inside the readout mode and x the dispersive coupling constant. The resonator 
frequency is thus shifted by ±x/27r = ±1.75[— 0.11/ + 0.14] MHz depending on the TLS 
state (see supplementary informatioijB| . The phase (f oi a microwave signal at Uc therefore 
acquires a TLS state-dependent shift after being reflected by the resonator, and provides a 



non-destructive readout of the TLS as demonstrated in numerous experiments 
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In our setup, the reflected signal is routed through a circulator to a cryogenic amplifier 
and is then measured by homodyne detection at room temperature, yielding the two field 
quadratures I{t) and Q{t). These time traces provide a continuous measurement of the 
TLS with a strength proportional to the signal input power, and thus to the intra-resonator 
average photon number n. In a fully quantum-mechanical description of the measurement 
process using the quantum trajectory formalism 19|], each quadrature can be written X{t) = 



X + {6X/2) {crz)c i't) +'Co(^)) where ((X^)^ (t) is the expectation value of a^ conditioned on the 
whole history of the detector outcome X{t') for t' < t, 6X is the maximum detector signal 
proportional to the measurement signal amplitude \/E, and ^o(^) is the total output noise 
of the amplifier. Our test of inequality (HI) consists in accurately measuring the steady state 
value of K{t) = {{X{t) — X){X{t + t) — X))t/{6X/2)'^ with a low measuring power, while 
the TLS is coherently driven. 

Measurement-induced-dephasing in ensemble averaged Rabi oscillations 

Recent experiments have already investigated the back-action of the measurement on a 
TLS with a similar circuit-QED setup [20]. However only ensemble averaged quantities (i.e. 
obtained by averaging the outcomes of many identical experimental sequences) had been 
considered prior to this work, and the only detectable effect of a measurement on the TLS 
dynamics was some extra dephasing, as demonstrated in 20| by measuring the broaden- 
ing of the TLS resonance line in presence of a field in the resonator. This measurement 
back-action results from the dependence of the TLS frequency Uge{t) = Uge + 2xn{t) on the 
photon number n{t) stored in the resonator: fiuctuations of n around n cause dephasing 
with a rate F^ (n) = Snx^ / k proportional to the measurement strength. We first perform a 
series of control measurements to verify on ensemble-averaged Rabi oscillations our quanti- 
tative understanding of the measurement-induced dephasing. After a field of n photons (see 
supplementary informatioifCl) is established inside the resonator using Vm-, a Rabi pulse of 
duration At is applied to the TLS with Vd (see Fig. [2j3), followed by a strong measurement 
pulse. The phase ((^(At)) of the reflected measurement pulse is averaged over an ensemble of 
typically 10"^ identical experimental sequences, yielding the data displayed in Fig. [2|3- One 
observes that for sufficiently low measurement strength n, the coherent dynamics is only 
weakly affected by the measurement. This is the regime where the Bell's inequality in time 



can be tested. For stronger measurement strength, the oscillations are progressively washed 
out and replaced by an exponential damping. For even stronger measurement strengths, 
the characteristic time of the exponential becomes longer and longer (see Fig. [2j3), revealing 
that a strong measurement inhibits the transition of the TLS from ground to excited state 
as expected from the quantum Zeno effect 2l|, |23|. We checked that these data are in quan- 



titative agreement with the expected measurement-induced dephasing (see supplementary 
informatioi j^ . However, it is important to realize that this set of measurements would be 
unchanged if our driven quantum TLS was replaced by a precessing classical spin, such as 
a macroscopic ferro magnet. Indeed, a macrospin obeys similar equations of motion as the 
expectation value of the spin of a TLS, namely Bloch equations. Thus no conclusion about 
the correlations between measurements at different times can be drawn that would allow a 
test of inequality (jlj) . 

Continuous measurement of Rabi oscillations in the frequency domain 

We measure these correlations by monitoring the system in its steady state, long after the 
transient ensemble averaged Rabi oscillations such as shown in Fig. 2B have been washed 
out. Instead of applying microwave pulses, the sources V^ and Vm are now continuously ON. 
The quantity of interest is the two-time correlation function K{t), whose direct calculation 
from the measured time traces X{t) is difficult in our setup because the amplifier noise 
dominates the output signal. However, this added noise can be removed by processing the 
signal in the frequency domain. 

For this purpose, we compute the square modulus Si and Sq of the Fourier transforms 
of I{t) and Q{t), to obtain the detector output power spectrum S{uj) = Sj{uj) + Sq{u). 
The signal power spectrum is then obtained by subtracting the amplifier noise spectrum 
Soff{^) measured when the two sources Vd and Vm are OFF from the signal-plus-noise 
spectrum Son{^) measured with both sources ON, and by dividing this difference by the 
independently measured frequency response R{u)) of the measuring line (see supplementary 
informatioi] [E|) . Typical curves are shown in Fig. [3l\. They show a single peak located 
at the Rabi frequency (already known from the time-domain measurements), without any 
harmonics within the 50 MHz detection window. The output spectrum of a continuously 
monitored TLS undergoing coherent oscillations 27H29| has been the subject of a number 
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Figure 3: Continuous monitoring of the TLS driven at the Rabi frequency w^j for different mea- 
surement strengths n. Each power spectrum is acquired in 40 to 80 minutes. (A) Spectral den- 
sities S'oAr(w) and Soff{^) when Vd and Vm are both OFF (brown) or both ON (red); here 
ujf(/2TT = 10 MHz and n = 1. The difference between ON and OFF shows a peak at the Rabi 
frequency. (B) Normahzed Rabi spectra Sz{u}) after correction from the frequency response of 
the measuring hne and conversion of the output voltage into units of cr^, at lo^J^tx = 5 MHz and 
n = 0.23, 0.78, 1.56, 3.9, 7.8, and 15.6 (from blue to red). The curves show the weak to strong 
measurement transition. (C-F) Normalized Rabi peaks at ujR/27r = 2.5, 5, 10 and 20 MHz (from 
blue to red) for n = 0.23 (C), 1.56 (D), 3.9 (E), and 15.6 (F). Thick and thin color lines are re- 
spectively the experimental spectra and those calculated from a theoretical analytical formula (see 
text and supplementary informatiorTG]) using only independently measured parameters (including 
x/27T = 1.8 MHz). Dashed black lines on top of the orange curves in (C,D,E) are Rabi peaks ob- 
tained by numerical simulation with the same parameters. The dotted-dashed black curve in C is 
the Lorentzian frequency response C(a;) of the resonator. 



of theoretical calculations 



i4, l25| , and precise knowledge of all sample parameters allows 
us to obtain for the first time a quantitative comparison with these theories. Indeed, we can 
convert the signal power spectrum in spin units by dividing it by a conversion factor {5V/2Y 
measured in a calibration experiment by saturating the g — e transition (see supplementary 
informatioi jFj) . The variation of the resulting Sz{<^) spectrum with increasing measurement 



power is shown in Fig. [3j3 and is in good agreement with theoretical predictions |9|, \2% 125] . 
These data clearly show the transition from weak to strong measurement in a continuously 
monitored driven TLS: at low n, the spectrum consists of a single Lorentzian peak at ur, 
upon increasing the measurement strength, the Lorentzian broadens towards low frequencies, 
and for strong measurements, the spectrum becomes a Lorentzian centered at zero frequency, 
similar to that of an incoherent TLS jumping stochastically between its two states. In 
terms of quantum trajectories, the Lorentzian spectra obtained in these two regimes are 
indirect signatures of the weak measurement-induced quantum phase diffusion along the 
Rabi trajectory, and of the quantum jumps made by the spin during the strong measurement. 
The theoretical curves shown in Fig. [HI are obtained using an analytical formula derived 
from the solution of Bloch equations [26] in which the finite detector bandwidth is taken 
into account phenomenologically (see supplementary informatioi fGj) ; the accuracy of this 
formula was checked by direct numerical integration of the system's master equation (see 
Fig. [3p-E anc fTTj) . The agreement between theory and experiment is good for n < 5 but is 
only qualitative at larger n, possibly due to a breakdown of the dispersive approximation. 

Experimental test of the Bell's inequality in time 

We now turn to the test of inequality ^. We measure a Rabi peak at bJ^jl-R = 10.6 MHz 
with n = 0.78 photons and a 30 MHz detection window. Under macrorealistic assumptions, 
the only effect of the bandwidth of the resonator would be to reduce the measured signal by 
its Lorentzian response function C{uj) = 1/[1 + {2uj/k,)'^]; we thus have to correct for this effect 
by dividing the measured spectral density Sz{uj) by C{u)). We then compute K{t) by inverse 
Fourier transform of Sz{u}) = S z{uj) / C (u) . The experimental and theoretical Rabi peaks as 
well as the corresponding /lg(t) curves are plotted in Fig.[4l showing good overall agreement 
despite residual low-frequency noise possibly originating from low-frequency fluctuations of 
Uge- The error bars on /lg(t) are the sum of the systematic errors in the calibration of 
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SV, K, and R{u)), and of the statistical error on the measured spectrum (see supplementary 
informationlj]). 
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Figure 4: Experimental violation of the "Bell's inequality in time" introduced in Fig. [TJ3. A) ex- 
perimental (red) and theoretical (blue) spectral densities 5*2, calculated or measured at uj /2Tr = 
10.6 MHz and n = 0.78. The experimental curve is obtained by correcting the raw 5*2 spectrum ( 
black line) , acquired in 13 hours with a 30 MHz bandwidth, from the frequency response C{u}) of 
the resonator (see Fig. [3p). The blue curve is calculated with F^ = 200 ns and F2^ = 150 ns (see 
supplementary informatior TGl) . B) experimental (dots) and theoretical (blue line) Leggett-Garg 
quantity fhciT) = '2K{t) — K{2t), with K(t) the signal autocorrelation function obtained by 
inverse Fourier transform of the Sz curves in the left panel. Green error bars correspond to the 
maximum systematic error associated with calibration and C{uj), whereas red ones also include a 
two standard deviation wide statistical error ib2(T(r) associated with the experimental noise on Sz- 
The Leggett-Garg inequality is violated (yellow region) at r = 17 ns (see green arrow) by 5a. 

We first note that we find K{0) = /lg(0) = 1.01 ± 0.15, a value close to 1 that directly 
results from the independent calibration of 6V. Since K{0) represents the variance of z{t) 
and \z{t)\ < 1, this confirms that at any time z{t) = ±1 as expected for a quantum TLS. 
A classical macrospin oscillating as shown in Fig. |3] and calibrated with the same method 
would have given a variance of 1/2 instead. Note also that K{0) = 1 could never be deduced 
from ensemble averaged Rabi oscillations such as those of Fig. |2|3- This is an example of the 
specific interest of correlation-function measurements compared to time-domain ensemble 
averaged signals. Most importantly, we observe that fhci'T') goes above the classical limit 
of 1, reaching 1.44(±0.12) ± 2 x (a = 0.065) at r = 17ns ~ tt/Sw/j and thus violating 
inequality (jl]) by 5 standard deviations a. This maximum of fici slightly below the ideal 
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value of 1.5 in absence of decoherence, is a direct signature of the invasive character of 
the measurement process, which projects partially but continuously the TLS towards the 
state corresponding to the detector output. It is the interplay between this continuous 
projection and the coherent dynamics that yields the violation of the inequality. More 
quantitatively, the maximum of f^c is in agreement with the quantum prediction of 1.36 
when taking into account the independently measured relaxation and dephasing rates of the 
TLS. This violation of the Leggett-Garg inequality rules out a simple interpretation of K{t) 
as the correlation function of a classical macrospin. It therefore brings further evidence 
that a collective degree of freedom characterizing a Josephson circuit can behave quantum- 
mechanically. It also demonstrates that the back-action of a weak measurement, far from 
being a simple noise that spoils quantum coherence as could be deduced from ensemble 
averaged measurements, tends also to reinforce correlations between measurements made at 
different times. 

It is interesting to discuss in what respect the assumptions made in analysing the ex- 
perimental data influence the final result: apart from simple corrections relying on classical 
electromagnetism, we determine the main normalization factor 6V by saturating the TLS 
transition and assuming that the ensemble averaged spin, either classical or quantum, obeys 
Bloch equations. We checked this assumption in Fig. [2] and supplementary Fig. [TJ which 
show in particular that the excursion of the signal when driving the spin is symmetric around 
the saturation value, as it would be for classical macrospins. The observed violation is thus 
not an artefact of our analysis framework, and represents more than a mere self-consistency 
check of a quantum model. 

Conclusion 

In conclusion, we have reported the experimental violation of a "Bell's inequality in time" 
by continuously monitoring the state of a superconducting artificial TLS while it performed 
Rabi oscillations. The measured two-time correlation function of the detector output reveals 
strong non-classical correlations between the signal already recorded and the TLS subsequent 
evolution. Our work thus brings a further proof of the truly quantum-mechanical character 
of Josephson artificial atoms. It is moreover a first step towards the test of a number of 
important predictions for such a system (9J, and towards certain quantum feedback schemes: 
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if the continuous monitoring could be performed with a quantum-hmited amphfier |3d.l3l|. it 
would become indeed possible to stabilize the phase of the Rabi oscillations by feeding back 
the demodulated signal onto the amplitude or frequency of the source that drives the TLS. 
This would modify the shape of the Rabi peak, on top of which a narrow line should develop 
[lO|. The correlations demonstrated in this work could then constitute a key resource for 
quantum feedback, in the same way as entanglement is a resource for quantum information 
processing. 
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Supplementary Information 

A. Sample fabrication, setup, and measurement protocol 

The sample is fabricated using standard lithography techniques. In a first step, a 200 nm 
Niobium thin-film is sputtered on a high-resistivity oxidized Silicon substrate. It is patterned 
with optical lithography followed by reactive ion etching of the Niobium to form the Coplanar 
Waveguide resonator. The transmon is then patterned by e-beam lithography followed by 
double-angle evaporation of two Aluminum thin-films, the first one being oxidized to form 
the junction oxide. The sample is glued on a microwave printed-circuit board, enclosed in 
a copper box, and thermally anchored to the mixing chamber of a dilution refridgerator at 
typically 20mK. 

Measurement signals are generated by mixing the output of a microwave source (V^ or 
Vm) with DC pulses generated by arbitrary waveform generators, using DC coupled mixers 
(not shown in Fig. 1). They are then sent to the input microwave line that includes bandpass 
filters and attenuators at various temperatures (77 dB in total). The output line contains 
a 4 — 8 GHz bandpass filter, a circulator and two isolators (not shown in Fig. 1), and a 
cryogenic amplifier CITCRYOl-12 (from Caltech) with 38 dB gain and noise temperature 
T/v = 4K. The output signals are further amplified at room-temperature yielding a total 
gain of 56 dB, and finally mixed down using an I/Q mixer with a synchronized local oscillator 
at the same frequency. The / and Q quadratures are filtered with a 50 MHz low-pass filter 
and further amplified with power gain 100 ; the total gains on channels / and Q are equated 
with a precision better than 0.5%. Both quandratures are then sampled by a fast digitizer 
and transferred to a computer that processes them. 

In the case of ensemble-averaged Rabi oscillations, 200 ns long measurement pulses are 
used; the corresponding X{t) traces are averaged over 10^ identical sequences (repetition 
rate: 200 kHz), yielding time traces as shown in Fig. [2]B. 

In the case of power spectrum measurements, we compute the Fast Fourier Transform 
X{uj) on records of 1024 X(t) samples separated by 10ns. We average Sx = \X{u)\'^ over 
10^-10^ identical sequences, yielding spectra as shown in Fig. 3A. In these experiments, it 
is crucial to accurately subtract the amplifier noise, whose power is about 60 times stronger 
than the Rabi peak we want to measure. For that purpose, we measure separately and 
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Figure 5: Experimental procedure for measuring the signal power spectra. Both microwave sources 
Vd and Vm are switched ON and OFF during Tqn = Tqff = 2.5 ms. Both quadratures are 
sampled, their Fast Fourier Transform calculated and squared to compute the power spectra Sx ON 
and SxoFF- Dashed rectangles show the Tss = 5;us waiting times for establishing the system 
steady state. 

subtract the signal+noise Sxon and the noise Sxoff by alternating periods of duration 
Ton = Tqff = 2.5 ms, during which both microwave sources V^ and Vm are ON, or both 
are OFF. Each time the pulses are switched ON or OFF, we wait a time Tgs = 5/is to let 
the system reach its new steady state. The resulting experimental sequence is shown in 
supplementary Fig. [5l Note that this subtraction also suppresses the influence of long-term 
drifts of the amplifiers gains. 

B. Sample parameters determination 

The resonator parameters Uc and k, given in the text, as well as the TLS frequency Uge 
are determined by standard spectroscopic measurements and are consistent with the design 
values. In all the experiments with no field in the resonator, the TLS frequency is fixed at 
Uge/27T = 5.304 GHz by tuning the magnetic field. The x value, which is used to calibrate 
the mean photon number n, is determined as follows: 

When the TLS state changes from g to e, the phase of the reflected microwave signal 
at the resonator frequency varies by 26ipo = 4arctan(2x//i:). We thus measure 6ipo by 
applying a strong and long driving pulse Vd that saturates the g-e transition and results in 
an almost equal 50% population of g and e. In this experiment, the resonator is probed 
continuously with a low-amplitude field (n ~ 2) while the reflected signal is measured by 
heterodyne detection with a local oscillator detuned from Vm by 3.2MHz. The resulting 
beating pattern in the two quadratures X{t) is ensemble averaged over a few 10^ identical 
sequences, and its phase (f{t) is fitted and plotted in supplementary Fig. El The plot yields 
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Figure 6: Ensemble averaged measurement of ip{t) when a driving pulse saturates the g — e transition 
starting from g. The measurement is continuous with n ~ 2 photons. 



a shift 6(j)Q = 12.5 ± 0.2° between no driving and saturation. 

For converting this shift into a x value, we take into account two effects: 

(i) the TLS is not perfectly in its ground state g in the absence of microwave drive due 

to residual thermal excitation to e. We have measured this thermal population p(e)o = 

0.02 ± 0.01 by measuring the noise spectrum with Vd being switched OFF. 

(ii) the saturation induced by the driving field is slightly below 50% due to longitudinal 

and transversal relaxation. The population p{e)st in the excited state at saturation is indeed 

given by the steady state solution to Bloch equations 



p{e] 



1 .1 



St 



p{e) 



1 + {T^'6uy 



1 + (r^-^^c.)^ + {ulT^^T 



2r-ir-iA' 



(5) 



^-1 

2 



Here ujr/2tt = 10 MHz is fixed by the driving strength, Fj"^ = 200 ± 10 ns and F2 
150 ± 10 ns are independantly measured, and 6u) is the residual detuning of the driving 
source from the TLS resonance. In our experiment the driving frequency was scanned to 
experimentally maximize 50o so that 6u ^ . This yields p{e)st = 0.496 ± 0.001 instead of 
0.5. 

Taking (i) and (ii) into account, we thus calculate 6(f)o = (0.95 ± 0.02)6ipQ. Maximizing 
all the uncertainties including those on k, we finally obtain x/27r = fi;tan((5(y9o/2)/47r = 
1.75(— 0.11/ + 0.14) MHz, also consistent with the TLS parameters determined by spec- 
troscopy. 
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C. Calibration of the mean photon number n at frequency ooc 

The cahbration of rl as a function of the apphed input power P of Vm is needed for 
checking that the measurement-induced dephasing is quantitatively understood (see D) and 
for producing the theoretical curves of Figs. |3] and H] without any fitting parameters. To 
perform this calibration, we measure the AC Stark shift— 2x?T of the TLS frequency uoge 
as a function of the input power P; x being known (se jBJ) . the shift provides an in-situ 
calibration of n{P). 

Note that all the curves presented in this article were obtained with the TLS driving 
source Vd tuned in resonance with the AC Stark shifted TLS frequency. For that, the TLS 
spectroscopy peak in presence of the very same n field was always measured just before each 
acquisition. 

D. Measurement-induced dephasing and Quantum Zeno Effect 

Since measurement-induced dephasing is a key ingredient of the present work, we checked 
our quantitative understanding of it. For that, we measure ensemble-averaged Rabi oscil- 
lations (see supplementary Fig. [7]) in presence of a perturbing field of n photons at fre- 
quency Uc, which mimicks the effect of a continuous measuring field. In this experiment, 
the AC stark shifted TLS frequency ujge{t) = ^ge — '^X^i't) fluctuates with the number of 
photons n{t), which fluctuates because of the shot noise in the measurement signal. In an 
ensemble-averaged viewpoint, this leads to pure dephasing of the TLS with a dephasing rate 



r^ (n) = Snx^ / K already observed in other experiments 20|. Now, because Rabi oscillations 
at ujr are most sensitive to the noise spectral density at ur S3[ and because the shot noise 
in the measurement signal is flltered by the resonator response C{uj) = 1/[1 + {2u/k)'^], one 
expects for the Rabi oscillations a pure dephasing rate 

Tf{u:n,n) = Tf{n)C{u:^). (6) 

Supplementary Fig. [7| conflrms all these expectations: the Rabi oscillations are progres- 
sively washed out and replaced by an exponential damping when the measurement strength 
n is increased. For even stronger measurements (larger n), the damping time constant in- 
creases (see Fig. EB), indicating an inhibition of the TLS transition from g to e that is a 
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Figure 7: Determination of the measurement-induced dephasing using ensemble-averaged Rabi 
oscillations recorded with the protocol of Fig.[2j3. (A) Same data as in Fig. [2j3 with uj{^ /2ti = 5 MHz 
and n = 0, 1, 2, 5, 10 and 20 (blue to red). Fits to the solution of Bloch equations (thin black 
lines) yield the measurement-induced dephasing rate F^ , after subtraction of other decoherence 
contributions deduced from the curve at n = 0. (B) Similar experimental data at uJii/2'K = 2.5 MHz 
and same photon numbers. When the measurement strength is increased, the time dependence at 
short times changes from quadratic to approximately linear with an increasing time constant, which 
is a manifestation of the Quantum Zeno effect. (C) Experimental values (dots with error bars) of 
F^ (n) for ijjfi/2TT = 2.5, 5, 10 and 20 MHz (blue to red). Solid lines are linear fits of the data taken 
at n < 5 , yielding the measurement-induced dephasing rate per photon 7^ . (D) 7^ as a function 
oiuju. Comparison between the experimental values (magenta dots), values obtained by numerical 
integration of the system master equation (blue squares), and theoretical curves (magenta lines) 
given by Eq. [6] and using only measured parameters (the two lines limiting the grey area correspond 
to the lower and higher bounds of experimental uncertainties) . 

signature of the Quantum Zeno Effect. By fitting each Rabi curve with the analytical solu- 
tion of Bloch equations 26| (using the measured F^^ = 225 ± 10 ns), we obtain in Fig. [7p 



the total decoherence rates T2{ijJRin) that includes the measurement-induced dephasing 
F^ {ujR,n) and the contribution T2{cor,0) = T^{ujfi) + Fi/2 from other dephasing sources 
and from energy relaxation (F°(27r x 5 MHz) = (810±20ns)~^ in Fig. [7JA. for instance). For 
each cor, the measured F^ {cor, n) are proportional to n as expected. Their slopes 7^ {^r) 
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determined by fitting F^ (cor, n) up to n = 5 are shown in Fig. [7p and are in good agreement 
with the predictions of Eq. [6] using the independendly measured values of x and k. We thus 
have a quantitative understanding of the measurement-induced dephasing in our system. 

E. Determination of the frequency response R{uj) of the measuring line 



The measuring hne between the sample and the output signal includes several microwave 
circulators, filters, amplifiers, an IQ demodulator, two fast digitizers, and many sections 
of cables with connectors. Its frequency response R{u) was measured in-situ in order to 
correct precisely the raw measured spectra S{uj). It is shown in supplementary Fig. [S] after 
normalization to 1 at zero frequency. The error bar shown (±1.5%) is an upper bound of 
the maximum systematic error over the whole frequency range. 
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Figure 8: Frequency response R{u)) of the measuring line, including the amplification and demod- 
ulation chain. The error bar represent a constant maximum relative error. 



F. Calibration of the spectra : dV{n) 

According to the definition of S{u)), the conversion factor 6V{n)/2 between the demod- 
ulated output signals in Volts and spin units is defined as [(5V(n)]^ = [5/(n)]^ + [SQln)]"^, 
with 6X(n) the change in quadrature X when the TLS state changes from g to e. This 
definition has the great advantage of being insensitive to any drift or jitter of the rela- 
tive microwave phase between the measurement source Vm and the local oscillator used 
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Figure 9: Calibration of the conversion factor bV l2 between output demodulated voltage and spin 
units. A measurement pulse Vm is applied with an amplitude corresponding to n = 0.78. After 
about one jis, a saturating pulse is applied at w^e during about 2 ps. The V'^ii) signal is measured 
and averaged over a few 10^ identical sequences. Starting from a thermal mixture of 98% g and 
2% e, the TLS undergoes Rabi oscillations at about 20 MHz before it reaches its steady state with 
50% g and e population. The steady-state output yield SV (see text). 

for the demodulation. Since [bV jl)^ is the normalization factor of the spectrum used to 
test inequality (4), it is particularly important to determine it with the best possible pre- 
cision. Although it could be calculated from several independently measured parameters, 
we found more accurate to calibrate it by direct measurement: we thus ensemble average 
V^if) = [IoN(t) — loFFit)]"^ + IQoNit) — QoFFit)]"^ Under saturation of the g-e transition 
(see supplementary Fig. [9]), similarly to what is reported in B to determine x- In this ex- 
periment, Vm is always ON and only Vd is switched ON and OFF, and ujr/2'k = 20 MHz 
so that p{e)st = 0.499 ± 0.001; V^{t) varies from e = 4 < .^q > (with < ^q > the vari- 
ance of the noise on each X up to f ^ + e with v = [p{e)st — p{e)o]SV when saturation is 
reached. From supplementary Fig. [9l we obtain {6V)'^ = 10.29 ± 0.64 mV^ for n = 0.78. It 
is important to note here that this calibration was done with the very same measurement 
power as that used to record the spectrum of Fig. |U so that the uncertainty on n does not 
impact at all the violation of the Leggett-Garg inequality. Moreover, as averaging of this 
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particular spectrum took about 13 hours, the cahbration was performed twice, before and 
after averaging, in order to check that the experiment was stable: the second calibration 
yields {6Vy = 10.59 ± 0.64 mV^. We thus take {6V/2y = 2.61 ± 0.16 mV^ for the conversion 
factor in our test of inequality (jl]). 

The experimental spectra of Fig. 3 obtained for different n have also been expressed in 
spin units by simply rescaling SV{n = 0.78) by y^n/0.78 at low n. For n above 3, we 
also take into account corrections to the dispersive approximation, which is valid only for 

n < n^rit = ^^l^gl = 31 
Using a model similar to 



a 



in our case (where go is the TLS-resonator coupling constant). 
19|, we keep the form of the dispersive Hamiltonian unchanged 
but use a modified dispersive constant x(^) = x(0)(l " ■^^) yielding a modified conversion 
factor 6V{n){l — An). We determine A = 7 • 10^^ both experimentally and theoretically. 
Note that this correction gives noticeable effect only for ra = 15 in Fig. [Sf. 

G. Analytical formula for the frequency spectra 

The noise spectrum of a TLS undergoing coherent oscillations under continuous mea- 
surement has been computed in the case of a quantum dot coupled to a quantum point 
contact. Such a problem can be completely mapped onto our experiment with a transmon 
in a resonator with infinite bandwidth. In this limit, an exact analytical formula exists 
for the spectrum 5*2(0;) [9], which is simply the Fourier transform of the two-time correla- 
tion function K^ij) = {cfz{t)(fz{t + T))t obtained from the analytical solutions to the Bloch 
equations: 






+ [(1 - 4)(r2 - rO/2 - colzl/T,] (f + Ql-u')}, (7) 



where ur = ^00% - (Fs - T,y/4, i = {T^ + T,)/2, and z,t = -1/(1 + T^'T^Wj^) is the 
steady state solution. Due to the finite bandwidth of the detector, the measured signal Sz 
is reduced as well as the dephasing rate according to Eq. (jH]). We modify Eq. [7] to take this 
into account by (i) multiplying Sz by the Lorentzian cutoff C{u) = 1/(1 + {2u/k)'^) and (ii) 
by changing the dephasing contribution F,^ to the total decoherence rate F2 = T^ + ri/2 into 
T^C{u); we obtain in this way an expression Sz{uj). The vahdity of this phenomenological 
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approach was checked by numerical simulations as explained in [Hi The calculated spectra 
of Fig. |3p-F and in the inset of Fig. [Hare obtained using the independently measured values 
Fi = (200ns)-iand T^C{uj) = (150 ns)-^ (see[D]). 



H. Numerical simulations 
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Figure 10: Comparison between the analytical spectra Sz (solid lines) and numerical simulations 
using the master equation as explained in the text (dashed black lines), for n = 0.23 and 3.9, and 
uptJ^TT = 2.5, 5, 10, 20 MHz. The agreeement is excellent. 

We perform numerical calculations of the combined TLS-cavity field evolution by inte- 
grating the master equation of the system (as described in [19]). We model its coherent 
evolution with the effective dispersive Hamiltonian (valid in the dispersive limit |A| ^ go) 



— ^^- Oz + {uc - Um)a)a' + x<^'^ci(^z + Vm{a^ + a) + Vd((T_ + o-+) (8) 



with ujd and uj^ the frequencies of the sources Vd and Vm driving the TLS and the resonator. 
Taking into account the damping of the resonator field at rate k and of the TLS energy and 
coherence at rates Fi and F^, the master equation is 



p = --[H,ff,pit)] + t,V[a]p{t) + Fil?[a_]p(t) + F^2)[a,]p(t)/2, (9) 

with T^IAIp = ApA^ — A^ Ap/2 — pA'^ A/2. We integrated this equation using a quantum 



optics library available online 



Sl| 
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We used this numerical tool to perform simulations of Rabi oscillations dephased by n 
photons in the cavity (dark blue squares in supplementary Fig. [Tli), taking x/2vr = 1.8 MHz. 
We also used it to compute the two-time correlation function of the detector output, i.e. 



K'{t) = K,{Re[a^{t + T)a{t)])t, by using the quantum regression theorem S2|]; we then 
convert K'{t) into a power spectrum by Fast Fourier Transform. The conversion into spin 
units was done as in the experiment, by calculating the output signal at saturation 6V of a 
continuously monitored Rabi oscillation. We were able to perform the simulations only up 
to n = 7.8, yielding the dashed lines shown in Fig. [Sp-F. The agreement with the analytical 
formula is excellent for all calculated curves, as can be seen in Fig. [TUl 

I. Experimental uncertainties in the determination of /lgI^) 

As mentionned in the text, the experimental points /lg(t) = 2K{t) — K{2t) of Fig. H] 
are obtained from the inverse Fourier transform K{t) of Sz{u)) = Sz{u))/C{u)) with 

oAUJ) — 



R{io){6V/2y ' 

with -R(O) = C(0) = 1 and 6V/2 being measured at zero frequency. The systematic error 
bars on fici'^) (in green in Fig. H]) result from the sum of the three maximum relative uncer- 
tainties, AR/R = ±1.5% , A{6V/2)y{6V/2f = ±6.1% and AC/C = 2(A/«//s:)/[l±(fi;/2a;)2] 
with Ak/k, = ±2.6%. Note that the frequency dependent error AC/C is propagated exactly 
through the calculation of /lg('^) and contributes to AfLc/fiG by ±0.8% where inequation 
dl]) is violated. The total systematic error at that point is thus ±8.4%. The statistical 
standard deviation on each fhci'T) data point is computed by propagating the statistical 
error on the measured Rabi spectrum (i) through the division by the cavity filtering C(a;), 
(ii) through the definition of the inverse Fourier transform, and (iii) through the difference 
2K{t) — K{2t). Each point A; = 1 to A^ of the Sz[iS) spectrum with bin size A/ = 100 kHz 
has a constant standard deviation a^ measured in the 22 — 30 MHz region where the spec- 
tral density is zero. Consequently, the standard deviation on each point k of the corrected 
spectrum is a^ = cro/C[2TTAf{k — 1)]. Finally, the standard deviation cr^ on each point r of 
heir = (r - l)/{NAf)] is 
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A conventional 2o"^, statistical contribution is added to the systematic error to form the total 
red error bars of Fig. HI At the second point r = llns where inequation (j31) is violated, 
the standard deviation s = (T2 = 0.065, and the bottom of the systematic error bar is 4.9 s 
above 1. 
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